Hydrodynamic model for electron-hole plasma in graphene 
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We propose a hydrodynamic model describing steady-state and dynamic electron and hole trans- 
port properties of graphene structures which accounts for the features of the electron and hole 
spectra. It is intended for electron-hole plasma in graphene characterized by high rate of inter- 
carrier scattering compared to external scattering (on phonons and impurities), i.e., for intrinsic or 
optically pumped (bipolar plasma), and gated graphene (virtually monopolar plasma). We demon- 
strate that the effect of strong interaction of electrons and holes on their transport can be treated as 
a viscous friction between the electron and hole components. We apply the developed model for the 
calculations of the graphene dc conductivity, in particular, the effect of mutual drag of electrons and 
holes is described. The spectra and damping of collective excitations in graphene in the bipolar and 
monopolar limits are found. It is shown that at high gate voltages and, hence, at high electron and 
low hole densities (or vice- versa), the excitations are associated with the self-consistent electric field 
and the hydrodynamic pressure (plasma waves). In intrinsic and optically pumped graphene, the 
waves constitute quasineutral perturbations of the electron and hole densities (electron-hole sound 
waves) with the velocity being dependent only on the fundamental graphene constants. 



I. INTRODUCTION 

The hydrodynamic approach is quite reasonable for the 
description of dense electron-hole plasma in semiconduc- 
tor systems in which the electron-electron, electron- hole, 
and hole-hole collisions dominate over the collisions of 
electrons and holes with disorder. In particular, such 
a situation can occur in intrinsic graphene at the room 
temperature, particularly, in the structures with low or 
moderate permittivity of the layers between which the 
graphene layer is clad, when the characteristic inter- 
carrier collision frequency can reach high values. As a 
result, this frequency can be much greater than that of 
collisions with impurities or phonons. Similar situation 
occurs in the gated graphene at sufficiently high gate 
voltages when large electron or hole densities can be in- 
duced as well as at strong optical pumping of graphene. 
The hydrodynamic models of the electron-hole systems 
in different structures and devices based on the standard 
semiconductors with parabolic and near-parabolic energy 
spectra of electrons and holes are widely used (see, for in- 
stance, Refji). However, in case of graphene, such mod- 
els and the pertinent equations should be revised due to 
the linear energy spectra of electrons and holes. A hy- 
drodynamic approach was recently used to describe the 
stationary transport processes in graphene'^. However, 
the role of a strong electron-hole scattering, which, as 
shown in our work, is crucial, was not addressed. 

In this paper, we develop a strict hydrodynamic model 
for electron-hole plasma in graphene and demonstrate 
its workability in some applications. The paper is or- 
ganized as follows. In Sec. II, we derive the hydrody- 
namic equations (continuity equations, Euler equations, 
and energy transfer equations) for graphene from the 



Boltzmann-Vlasov kinetic equations for massless elec- 
trons and holes assuming high carrier-carrier collision 
frequencies. Section HI deals with the application of the 
derived equations for the calculations of graphene dc con- 
ductivity taking into account the effect of electron-hole 
drag. In Sec. IV, the collective excitations in graphene 
at different conditions are considered using the obtained 
hydrodynamic equations. In particular, we show that 
two types of weakly damping excitations can exist in the 
electron- hole plasma: (1) electron (or hole) plasma waves 
in gated graphene in the state with the Fermi level far 
from the Dirac point and (2) quasi-neutral electron-hole 
sound waves in the bipolar electron-hole plasma. The 
damping of both plasma waves (in virtually monopo- 
lar plasma) and electron-hole sound waves (in bipolar 
plasma) is determined by the scattering on disorder, 
while the electron-hole scattering almost does not affect 
the damping. In contrast, the plasma waves in bipolar 
electron-hole systems exhibit a strong damping due to 
electron-hole scattering processes. In Sec. V, we draw 
the main conclusions. Some cumbersome formulas are 
singled out to the Appendix. 



II. DERIVATION OF HYDRODYNAMIC 
EQUATIONS 



For the massless electrons in holes in graphene the 
spectrum is linear e{p) = vpp. Hence, the kinetic equa- 
tions governing the distribution functions /e = /e(p) and 
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fh = fh{p) read, respectively, 



dt p dr dr dp 



St{f,J,} + St{f,,h} + SU{f,}, (1) 



dfh , pdfh 
at p or 



dr dp 

St{h, h} + St{h, fe} + SU{fh}- (2) 



Here vfp/p = de{p)/dp is the electron (hole) velocity, 
vp — 10^ cm/s is the characteristic velocity of electrons 
and holes in graphene (Fermi velocity), e = |e| is the 
absolute value of the electron charge, E = —dip/dr is 
the electric field, Sti{fe} and Sti{fii} are the collision 
integrals of electrons and holes, respectively, with dis- 
order (impurities and phonons); S'i{/e,/e}, St{fhJh}, 
and St{fi,, ffi} are the inter-carrier collision integrals. 
In Eqs. ([1]) and ([2]) we have neglected the recombina- 
tion terms assuming that the recombination rate is much 
smaller than that of collisions and the frequency of the 
plasma waves under consideration. 

As known^, the Fermi distribution functions of elec- 
trons and holes 
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turn the electron-electron and hole-hole collision integrals 
to zero owing to the conservation of momentum and en- 
ergy in the inter-carrier collisions. In Eqs. ^ and (|3]) 
Ve and V/i are the average (drift) velocities of electrons 
and holes, respectively, fie and iih sue electron and hole 
chemical potentials, and the temperature T is measured 
in energy units. 

Hereafter we consider rather small drift velocities and 
perform an expansion of Eqs. ([3]) and Q over Vg and V/ji 
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Here /g.o and fh.o stand for the functions, given by 
Eqs. © and (U), with Ve = V/j = 0. The distribution 
functions (j5l6p still turn the electron-electron and hole- 
hole collision integrals to zero and, at the same moment, 
conserve the number of particles (the density is the same 
in moving and stationary frames). 

If the electric field is sufficiently weak and the inho- 
mogenity of the electron-hole plasma is small (as it will 
be assumed in the following), the electron-hole collision 



integrals can be expanded over Vg and V/^ and presented 

as StifeJh] = -St{hJe} {Vh - Ve) • S{/e,o,A.o} 

(for details see the Appendix). In this form, the electron- 
hole collision terms describe the friction between the elec- 
tron and hole plasma components. Similarly, the terms 
corresponding to the collisions of electrons and holes 
with disorder can be presented in the form Sti{fe} — 
-Ve • S,{/e,o} and SUifh} = -Yh ■ S,{A,o}. Here S 
and Si are the functionals of the distribution functions 
/e,o and fhfi. 

At small values of electron and hole average velocities 
Ve and V/i, the friction terms can be considered as per- 
turbations in comparison with the electron-electron and 
hole-hole collision terms. Thus the distribution functions 
given by Eqs. ([3]) and ^ are the approximate solutions 
of Eqs. (HI) and Then the quantities Ve, V/i, /Xe, 

and iifi (or the electron and hole sheet densities, Ee and 
E/i) can be found considering the terms in the left-hand 
sides of Eqs. H]) and ^ and the friction terms as per- 
turbations, using the standard procedure (akin to the 
Chapman-Enskog method^ for the derivation of the hy- 
drodynamic equations from the kinetic equations). 

On integrating Eqs. ^ and ^ over dTp = 
g<Pp/ {2'iT%)'^ (where g = 4 is the electron degeneracy 
factor in graphene), we obtain the continuity equations 
for electrons and holes 
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To derive the Euler equations, one should integrate 
Eqs. H)) and ^ times p over dTp. In the case of 
parabolic dispersion, one could multiply the Boltzmann 
equation either by velocity or momentum as they are pro- 
portional to each other. For linear dispersion the choice 
of momentum is crucial. Just the momentum is con- 
served in particle-particle collisions unlike to the velocity 
which can be changed. After integration one obtains the 
system of Euler equations for electrons and holes: 

3 d (pe) Ve d Vp jPe) y _ 

2dt VF dr 2 ^ dr ~ 

- /3eVe - /3e/. (Ve - V^) , (8) 



3 d iph) V/i d vp iph) 



el^h — 



2dt VP dr 2 ' ^ " dr 

-/3hVh-fieh{Vh~Ve). (9) 

Here the angle brackets denote an integration over the 
equilibrium Fermi distribution functions, in particular. 
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is the momentum modulus per unit area, the friction co- 
efficients Peh: Pe and Ph are the functions of the non- 
perturbed (steady-state) values of the chemical potentials 
^e,o and fihfi- 
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One can rewrite Euler equations in a classical form on 
introducing the fictitious carrier masses 
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VF^h ' 



which are estimated as 0.016 of the free electron mass at 
= = and T = 300 K. 

The detailed derivation of the friction coefficients can 
be found in the Appendix; here we write down only the 
final expressions in several limits. The electron-hole fric- 
tion coefficient jS^h can be represented as 
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where k is the effective permittivity of environment (the 
substrate and gate dielectric), the dimensionless constant 
A of the order of unity and the function I {p^/T, fih/T) 
can be obtained after proper linearization of the electron- 
hole collision integral (see the Appendix). For intrin- 
sic graphene / (0, 0) = 1, while for monopolar plasma in 
gated graphene I {fi/T, ^/T) oc e^''^^ due to the expo- 
nentially small number of holes. In the current paper, 
we shall adopt the following interpolation for /3eh which 
is valid in the limiting cases of monopolar plasma and 
intrinsic graphene: 
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Here we have introduced the electron-hole collision fre- 
quency Veh to be estimated below in the Sec. III. 

When the acoustic phonon scattering dominates, the 
coefficients (3^ and fih can be presented as^ 



Ph = 



where D is the deformation potential constant, p is 
the sheet density of graphene, s is the sound veloc- 
ity. Due to the existence of several acoustic phonon 
branches and considerable discrepancy in experimental 
data of graphene constants^"— it is reasonable to use semi- 
phenomenological formulas: 



(12) 



and extract the numerical value of A from experimental 
data on dc conductivity (see Sec. III). 

For scattering on charged impurities, the friction co- 
efficients /3e and are reasonably proportional to the 
carrier densities and the density of charged impurities 
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In general, the density of charged impurities depends on 
the chemical potentials and temperature. 



At last, the system of hydrodynamic equations should 
be supplemented with the energy transfer equations. In- 
troducing the energy density (Sch) — VF{pe,h), the per- 
tinent equations can be written down as 
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where the heat sink rates Qe and Qh are proportional 
to the difference between electron and phonon temper- 
atures. However, in the present paper we neglect the 
heating as we deal with low-field dc conductivity and 
high-frequency plasma waves. 



III. EFFECT OF ELECTRON-HOLE DRAG AND 
DC CONDUCTIVITY OF GATED GRAPHENE 

As one of the demonstrations of the hydrodynamic 
equations applications we calculate the dc conductivity 
of gated graphene. It is assumed that the electric po- 
tential and the charge density are related by the Poisson 
equation. In the gradual channel approximation^ 



C{Vg -<p)^ e(Ee - S,0 



(16) 



Here C — K/And is the specific capacitance per unit area, 
d and k are the thickness and permittivity of the gate 
dielectric, respectively, and Vq is the gate voltage. 

To calculate the dc conductivity we rewrite the Eu- 
ler equations for the steady-state situation in terms of 
electrochemical potentials: 



d {fie - ef) 
dr 

djeip- Hh) 
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- -/?eVe - I3eh (V e - V,,) , (17) 



^PhVh - Peh {^h - Ve) , (18) 



where the derivative of electric potential is associated 
with drift current, as well as the derivative of chemical 
potential is associated with diffusion current. In the fol- 
lowing we restrict our calculations by the consideration 
of the drift current only and, therefore, omit the terms 
dfj,e,h/dr in Eqs. (|17|) and (fT8| . As a result, we arrive at 
the following expressions for the mean (drift) velocities 
of electrons and holes: 
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Substituting the quantities Vg and V/j from Eqs. ((12121) 
into the general expression for current density j = 



e{T,hVh — 5]eVe) = G {—d(p/dr), one obtains the fol- 
lowing formula for the conductivity G: 



G 



(21) 



The first term in the right-hand side of Eq. (pij) is 
associated with the scattering of electrons and holes on 
impurities and phonons. This term turns into zero at 
the Dirac point. However, far from the Dirac point, i.e., 
in purely electron or hole plasma, this term dominates. 
Meanwhile, the second term is due to the contribution 
of the electron-hole friction. The latter results in a high 
resistivity of graphene at the Dirac point and its vicinity. 

To calculate the graphene minimal conductivity, we 
reasonably assume that at the Dirac point scattering be- 
tween electrons and holes prevails: ^ Peh, Ph Peh- 
Plugging Eq. ([TU| for (3f.h into Eq. (PT|) one obtains an 
expression for intrinsic graphene conductivity 
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FIG. 1. Conductivity of graphene G vs. gate voltage at 
T = 300 K, d = 10 nm, and k = 4. Insets: left panel - 
resistivity G^^ near the Dirac point vs. gate voltage, right 
panel - mobility B vs. chemical potential. 
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Strikingly, the intrinsic graphene conductivity (minimum 
conductivity) does not depend on temperature. This is in 
agreement with the experimental results^, demonstrating 
a constant conductivity over a broad range of tempera- 
ture from 0.3 K to 300 K in which the carrier density 
varies by 6 orders of magnitude. Worth mentioning the 
Eq. (1^^ was previously obtained using the scaling the- 
orjiiS, and also via thorough description of electron-hole 
scatteringii. 

Earlier, the minimum conductivity of graphene was 
calculated under the assumption of strong interaction 
among carriersi^ril. The strength of interaction is gov- 
erned by the "fine-structure constant" a = e^/Khvp 
which is equal to 2.2 for intrinsic (suspended) graphene 
(k — 1). The strong-interaction theories bind the con- 
ductivity of graphene to the conductance quantum jh 
regardless of the permittivity of environment. As the 
parameter a is about unity, both theories of weak and 
strong interaction give rise to close values of conductivity. 
However, in prospective graphene structures the screen- 
ing caused by dielectrics and nearby gates make the value 
of a less than unity, therefore, the weak-interaction the- 
ories become more adequate to the situation. 

After substituting Eq. for Peh into Eq. (I^ one 
can rewrite that equation in the conventional form ap- 
propriate for the estimation of the electron-hole collision 
frequency v^h'- 
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where Mq is the fictitious mass of electrons and holes 
in intrinsic graphene dependent on temperature. Com- 
paring Go given by Eq. (j23|) with the experimental value 
Go — (6 kOhm)~^, one can estimate the electron- hole col- 
lision frequency as v^-h = 3 x 10^^^ . Such a high value 
justifies an employment of the hydrodynamic model. 



In the monopolar limit (Eg ^ E^), the expression for 
the conductivity could be also simplified: 



G = 



(24) 



Equation (|24l) allows to estimate the coefficient A in the 
phonon collision term given by Eq. (1121) . It is widely 
assumed that in suspended graphene samples or twisted 
graphene stacks the conductivity and mobility are limited 
by phonon scattering only. For instance^^, the carrier 
mobility B, defined as 



B = 



G 



reaches the value of 120,000 cm^ V"^ s^^ at Eg - E,, = 
2 X 10^^ cm~^ and the temperature T — 240 K (when the 
chemical potential /i ~ 0.028 eV). Hence, in the monopo- 
lar limit, one obtains 
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eE, 



(25) 



Equation ^ yields A ~ 2.4 x lO^'' cm s'^ J"2_ Ac- 
cordingly, the characteristic collision frequency z^e — 
2vp/3e/3{pe) at room temperature varies from 8.6 x 
10" s"^ at the Dirac point to 3.5 x lO^^ g-i Vb = 10 V 
and d = 10 nm. 

Figure [1] demonstrates the dependence of the low-field 
dc conductivity G of graphene on the gate voltage Vg, 
plotted using Eqs. ([2T1 [TTl [T2|) . The voltage dependence 
of the graphene resistivity and the dependence of the 
carrier mobility on the chemical potential /i [given by 
Eq. (pSj)] are shown in the insets in Fig. [T] In the case 
when scattering on acoustic phonons dominates, the con- 
ductivity tends to saturation at high gate voltages and 
the phonon-limited mobility decreases. To describe the 
conductivity of realistic graphene structures, one should 
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account for other scattering objects: charged impuri- 
ties, interfacial phonons, bulk phonons, etc. For exam- 
ple, scattering on charged impurities results in the fol- 
lowing trends of conductivity curves in monopolar case: 
G(Vg) oc Vg if the density of impurities is constant and 
G{Vg) oc \/Vg if the distribution of impurity levels is 
uniform in an energy scale^^. The latter fact manifests a 
decrease in defect scattering at low temperatures. 

The conductivity and mobility curves in Fig. 1 exhibit 
a good agreement with the experimental works, where 
phonon- limited mobility was investigated^^^-'^^. As for 
the theoretical approaches to the description of graphene 
conductivity, significant efforts have been focused on the 
calculation of the minimum conductivity and the conduc- 
tivity far from the neutrality point^iiii^. The conductiv- 
ity at an arbitrary value of carrier density was considered 
in Refs.^'*'^^ without specifying the mechanism of scatter- 
ing. Our model provides an opportunity to calculate the 
transport characteristics of graphene at any charge car- 
rier density taking into account the definite scattering 
processes. Above all, we show that the effect of electron- 
hole drag significantly impacts the transport properties 
of graphene and leads to an abrupt drop of graphene re- 
sistivity outside of the Dirac point. 

It is remarkable that the drag effect is crucial even at a 
small mismatch in the electron and hole densities^ . The 
conductivity in the regime of drag almost does not de- 
pend on electron-hole scattering, although it still remains 
the strongest one in the system. Indeed, the phonon (im- 
purity) scattering term in the expression for conductiv- 
ity (pij) dominates under the condition 

V P-i- ^=0 

For the preceding estimations, the ratio given by Eq. ([25)1 
is evaluated as 0.3 at room temperature. Providing 
ineq. (|26p is satisfied, the difference between electron and 
hole velocities (Eq. becomes rather small: 

\Vh-Ve\ ^ 1 . . 

» u— I 



The difference between velocities still exists because of 
the opposite charges of electrons and holes and, there- 
fore, the opposite directions of the electric forces. In 
principle, for the drag regime the two-fiuid description of 
graphene system might be reduced to a single electron- 
hole fluid with renormalized charge of particles according 
to Eq. (1^51) . However, it seems preferable to retain two- 
fiuid model to cover all possible situations. 

It is readily seen from above mentioned that a narrow 
resistivity peak in Fig. 1 could be explained by the drag 
effect. Moreover, stronger scattering leads to a wider 
resistivity peak. The latter is in concordance with the 
experimental data-'^''. 

IV. PLASMA AND ELECTRON-HOLE SOUND 
WAVES IN GRAPHENE 

A. General dispersion relation for collective 
excitations 

Below we demonstrate the application of the hydro- 
dynamic equations for the calculation of the spectra of 
collective excitations in electron-hole system. We use the 
same gradual channel approximation as in the previous 
section [Eq. p^ ]. To derive the spectra, one can apply 
the common technique of small-signal analysis assuming 
that 

Aie = /^e,o + Vee'('"""*^ , f^h = M/.,0 + S fihe'^''^~^'^ , 

where 6 fie, 6fih, '^V'e.o, ^Vh.o, and 6(p are the amplitudes 
of alternating variations. The non-perturbed chemical 
potentials fj-efl and fj,hfl coincide if the steady state is an 
equilibrium one and are determined by the doping and 
the gate potential. In particular, at Vg = and zero 
doping the Fermi level is located in the Dirac point. In 
intrinsic graphene under optical interband pumping the 
chemical potentials of electrons and holes can be rather 
different. For example, at Vg = 0, fj,e,o = Mo > and 
fih = —fJ-o < 0; where /iq is determined by the intensity 
of pumping. 

From the linearized versions of Eqs. ([7][9]) as well as 
Eq. (|16p the following system of algebraic equations 
arises: 
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The solvability condition for Eqs. (|28|) and ([29)1 results in the general dispersion relation for the collective excitations: 



-iuj + 



(1 + re) +Vf,+VF 



2f3eh 



-iuj H ^ (1 + rh) + Vh + VF 



iPe) . 
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Here for brevity we have introduced the dimensionless 
constants 



rh 



Cv p Cvp 
the squared characteristic velocities 



3 ipe) {Pe') ' 



and the frequencies 



3(Pe) 



^{Ph) {Pny 



HPh) ' 



The quantities Tg and determine the ratio of the elec- 
trostatic energy to the kinetic energy in the wave. In 
the particular gated structures we have considered above 
(k = 4, c? = 10 nm) the dimensionless constant varies 
from 1.9 at Vfe = to 45 at Vq = 10 V. The electron- 
phonon collision frequencies Ve and Vh, as shown below, 
determine the damping of the waves. 



B. Analytical solutions for symmetric bipolar and 
monopolar systems 

Exact solutions of the general dispersion equation ([50)1 
can be acquired for symmetric bipolar plasma (optically 
pumped or intrinsic graphene) and monopolar plasma. 

In symmetric bipolar plasma all quantities character- 
izing the electron and hole systems coincide (e.g. — 
Vh = f, Se = = E, ...), therefore, we shall omit 
the subscripts in this certain case. In symmetric systems 
Eq. ([50)1 provides two solutions: 
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(32) 

The second solution oj^ represents the waves with op- 
posite motion of electrons and holes leading to the strong 
damping. These waves correspond to the perturbations 
of charge density (plasma waves). 

The branch given by Eq. (PT|) is of the most inter- 
est. Electrons and holes in these waves move in the same 



direction, i.e., the electron- hole plasma remains quasi- 
neutral. Due to the co-directional motion of the electron 
and hole components, their mutual collisions do not af- 
fect the wave damping. In the case, it originates from 
the collisions of electrons and holes with impurities and 
phonons. We can conclude that those waves are asso- 
ciated solely with the pressure gradient. In the follow- 
ing, such waves will be referred to as electron-hole sound 
waves in analogy with electron-ion sound waves in clas- 
sical plasma. 

The velocity of electron-hole sound is 



2E2 



S_ = V = Vpi, 



3 ip) {p- 
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It does not depend on structure parameters and does 
not exceed the Fermi velocity for any electron and hole 
densities. If the Fermi level crosses the Dirac point the 
analytical expression for velocity is 



s_ = Vp 



18v/ln2C(3) 



0.6vp, 



where (^{x) stands for Riemann zeta function. 

The dispersion laws for monopolar plasma can be de- 
rived providing the inequalities Sg ^ ^h, Pe ^ Ph, and 
re ^ 1 ^ r/i are satisfied: 



(34) 



Veh 



V 2 3 



(35) 

A detailed analysis of these solutions shows that the 
branches aj+ and a;_ correspond to the oscillations of 
majority (electrons) and minority (holes) carriers, re- 
spectively. The minority carrier oscillations are strongly 
damped by electron-hole friction and, therefore, they are 
omitted in the further consideration. 

The wave velocity of the electronic plasma oscillations 
s-u is 



t;e(l + re) 



1/2 



(36) 



where a — e /nTivp is the coupling constant ("fine- 
structure constant") and kp = iJLf./hvp is the Fermi mo- 
mentum of electrons. This velocity markedly exceeds the 
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FIG. 2. Dispersion of plasma and electron-hole sound waves 
at difTerent gate voltages. Inset: non-linear dispersion of the 
waves at small frequencies resulting from collisions 



Fermi one and is primarily determined by the Coulomb 
interaction of carriers with the gate. 

The dispersions for electron-hole sound in intrinsic 
graphene and plasma waves in gated graphene at high 
gate voltages are depicted in Fig. 2 in the THz range of 
frequencies. The curves are calculated under the assump- 
tion of equal chemical potentials ^e,o = M/i,Oj i-e- when 
the non-perturbed state is an equilibrium one. The spec- 
tra exhibit non-linear dispersion at low frequencies orig- 
inating from the damping (see the inset in Fig. 2). 

Worth mentioning the dispersions are also non-linear 
at short wavelengths [kd ~ 1) when the gradual chan- 
nel approximation becomes inapplicable. However, this 
limitation can be overcome if we use the rigorous solu- 
tion of Poisson equation^! instead of Eq. ([16]). For the 
considered gated structure the transition to the rigorous 
solution can be performed by the substitution: 



C 



2Kk 



47r(l 



o-2kd 



(37) 



Equation (j37p restrains the unlimited growth of the 
plasma wave velocities for large distances d between 
graphene sheet and metal gate. 



C. Velocities and damping rates of the waves 

The wave velocities and damping rates can be calcu- 
lated analytically in the linear domain of spectra, i.e. at 
{i^e,h < w < s/d) at an arbitrary value of electron and 
hole densities. Assuming uj = sk, where s is the wave 
velocity, and plugging this into Eq. (pO)) , one obtains two 
solutions s_ and s+: 



1 



4 = 7rK(l + '^e)+<(l + r,0]± 



_|_ /^y^y'^Terh, 



(38) 



FIG. 3. Velocities of plasma waves vs. gate voltage calculated 
for different gate layer thicknesses [Eq. [36] . Dash-dotted line 
corresponds to the electron-hole sound velocity in the vicinity 
of the neutrality point. Regions of strong damping are filled 



consistent with Eqs. ([33l [36]) . 

In the vicinity of Dirac point the waves with the lower 
velocity s_ correspond to the electron-hole sound, while 
in the monopolar plasma they turn into the oscillations of 
minority carriers. The waves with the higher velocity s+ 
behave as plasma waves at any electron and hole densi- 
ties. The solutions given by Eq. ([55]) are plotted in Fig. 3 
as the functions of gate voltage with the assumption 
of equal steady-state chemical potentials (/ie.o = Mh,o)- 
One can readily see that the velocity s_ almost does 
not depend on the gate voltage applied, while the ve- 
locity s+ exhibits unlimited growth at high gate voltages 



1/4 



cx V. 



To calculate the damping rates of the waves we assume 
u)± = s±k -\- 17-1-, where 7± characterizes the damping 
rate of the waves. The obtained damping rates 7± vs. the 
gate voltage are demonstrated in Fig. [4] It is clearly seen 
that in symmetric bipolar plasma the electron-hole sound 
branch a;_ exhibits weak damping. In the monopolar 
case (Eg ^ Tih) the plasma wave, corresponding to the 
oscillations of the majority carriers, is weakly damped. 
The minimum damping rate for the waves considered is 
of the order of 5 x 10^^ s~^ and is determined by electron- 
phonon (or hole-phonon) collision frequency if the only 
scattering mechanism is acoustic phonon scattering. At 
high gate voltages the damping coefficient 7+ grows lin- 
early in accordance with the expression for characteristic 
collision frequency u,,. The maximum damping rate of 
the waves is of the order of electron-hole collision fre- 
quency Veh- In accordance with these calculations, the 
regions of strong damping in Fig. [3] are marked with fill- 
ing. 
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FIG. 4. Damping rates for the two branches of spectrum: lj^ 
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D. Comparison with other models 

Plasma waves in gated graphene were discussed in 
Refs3i^ using the kinetic approach under assumption of 
coUisionless transport. As expected, formally setting the 
collision frequency zero in Eq. ([M)) . we obtain the spec- 
trum of plasma waves, derived previously ^'^^ Meanwhile, 
in the kinetic approach the existence of the electron-hole 
sound could not be predicted as it was assumed that the 
wave velocity should overcome Fermi velocity. 

The dispersion curves of the plasma waves obtained 
via the hydrodynamic and kinetic theories are compared 
in Fig. 5 for the same parameters of the gated struc- 
ture (k = 4, d = 10 nm), gate voltage Vg = 10 V and 
three different collision frequencies. The only difference 
between the curves is in the non-linear part of spectra 
originating from collisions. 

Plasmons in graphene were also intensively studied 
within the analysis of the polarizability function n((7, lS) 
in random-phase approximation^^"^^. In particular, the 
square-root cj cx \fk and linear uj oi k dispersions were 
obtained for the non-gated and gated graphene, conse- 
quently. To obtain the plasma wave dispersion for the 
non-gated graphene in hydrodynamic model we tend d 
to infinity in Eq. ^\ and plug it into Eq. ([M]). Thus, 
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FIG. 5. Comparison of dispersions for plasma waves, calcu- 
lated using kinetic model?^ (disregarding collisions with im- 
purities and phonons, solid line) and hydrodynamic model for 
Vg = 10 V, d = 10 nm, k = 4, and three different collision 
frequencies. 



plasma waves in non-gated graphene at high electron den- 
sities exhibit the following dispersion: 



,2 j,2 



(H-4aM-(y)'-^y. 



(39) 



In the low frequency limit the above equation is simplified 

^j2akkF, (40) 

and the result quantitatively coincides with that pre- 
dicted in^i. What concerns the result obtained in^ for 
plasmons in gated graphene, it coincides both with that 
obtained within kinetic approach^^ and, hence, with our 
resuh p4)) . 

The transition from almost linear to square-root 
plasma wave dispersion at different gate dielectric thick- 
nesses is shown in Fig. [5] for the given value of electron 
density Eg = 5 x 10^^ cm~^. In the short-wave limit 
kd^ 1 (providing the applicability of the hydrodynamic 
model) the dispersion of plasma waves is linear again 
with the velocity depending only on the fundamental con- 
stants of graphene. Worth mentioning that the electron- 
hole sound waves are insensitive to the gate and exhibit 
almost linear dispersion at any gated structure parame- 
ters. 



V. DISCUSSION OF THE RESULTS 

The hydrodynamic model of electron and hole trans- 
port in graphene which takes into account the linearity of 
carrier spectra has been developed. Strong interactions 
between electrons and holes result in mutual frictional 
forces proportional to the mismatch in electron and hole 
drift velocities. The interaction of carriers with acoustic 
phonons and impurities is also governed by a friction term 
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in Euler equations. The estimated electron- hole collision 
rate is greater than that between carriers and phonons 
(impurities) at room temperature. 

The model has been applied to the solution of two chal- 
lenging problems: dc graphene conductivity and spec- 
tra and damping of the collective excitations in graphene 
structures. 

The hydrodynamic equations derived provide an op- 
portunity to calculate the conductivity of graphene sheet 
in a wide range of gate voltages including the Dirac 
point, where the influence of electron- hole collisions on 
the charge transport is dominant. However, when there 
is even a small difference (proportional to the ratio of 
external scattering rate to that of carrier-carrier) be- 
tween electron and hole densities the majority carriers 
drag the minority ones. In this case, the conductivity is 
governed by scattering on phonons and impurities. The 
effect of drag is particularly pronounced in high-purity 
(suspended) graphene samples. 

The spectra of collective excitations in electron-hole 
system in graphene have been calculated. The ex- 
istence of two types of the excitations has been re- 
vealed: quasineutral electron-hole sound waves and 
plasma waves. 

In the vicinity of the Dirac point or in the optically 
pumped intrinsic graphene the sound waves undergo 
weak damping of the order of 5 x 10"'^"'^ s~^ in perfect 
structures. The damping of such waves is caused by 
weak scattering on acoustic phonons and is insensitive 
to strong electron-hole scattering. The latter fact is due 
to co-directional motion of electrons and holes in quasi- 
neutral sound waves, which almost eliminates electron- 
hole friction. 

The quasi-neutral electron-hole sound waves under 
consideration are akin to those in bipolar electron-hole 
plasma in semimetals considered by Konstantinov and 
Perel a long time agc?^. Later such waves were stud- 
ied in semiconductors^^, two-band metals-^ and super- 
conducting compounds^S,. Recently the acoustic plasma 



waves were studied in two dimensional semimetalsSS, (Cd- 
HgTe/HgTe/CdHgTe quantum wells). It was shown that 
those waves can contribute to the thermodynamic prop- 
erties of the materials, e.g. carrier relaxation rates in 
bulk semiconductors^^ and critical temperatures of su- 
perconducting compounds^'^ . In graphene such waves 
could be feasibly excited by optical spots or modulated 
optical radiation. In both situations quasi-neutral non- 
uniform distributions of electrons and holes arise. 

The propagation of plasma waves in symmetric 
electron-hole systems is strongly suppressed by electron- 
hole collisions. Strong damping of such waves can sig- 
nificantly decrease the rate of electron-hole recombina- 
tion via plasmon emission, discussed in Refj^ In its turn, 
the recombination via the emission of electron- hole sound 
waves is prohibited as their velocity is smaller than vp 
[Eq. dSHl)]. 

In sufficiently monopolar systems (electron or hole 
plasma, e.g. in gated graphene at high gate voltages) the 
damping of plasma waves, corresponding to the oscilla- 
tions of majority carriers, is associated solely with the 
scattering on disorder (phonons and impurities). It can 
be rather weak in perfect structures, particularly at low 
temperatures. The dispersion of plasma waves in gated 
graphene structures is linear in a wide range of frequen- 
cies^S,. The dependence of the plasma wave velocity on 
the gated structure parameters is quantitatively similar 
to that obtained for the two-dimensional electron gas of 
massive electrons^. 

The hydrodynamic and kinetic models yield the same 
dispersions of plasma waves if one formally tends the col- 
lision frequency to zero. However, the hydrodynamic ap- 
proach provides a regular way to describe the damping 
of the waves, associated with carrier scattering. 

In conclusion, the hydrodynamic equations for bipolar 
graphene system were derived and applied for calculation 
of dc conductivity and spectra of collective excitations in 
graphene structures. The model opens the prospects to 
the simulation of graphene-based transistors, THz-range 
detectors and generators, and light emitting devices. 
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APPENDIX 1. DISSIPATIVE TERMS IN THE 
EULER EQUATIONS 

In this appendix we derive explicit expression for the 
friction forces, associated with carrier-carrier, carrier- 
impurity, and carrier-phonon scattering. 
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A. Mutual electron-hole friction 

To derive the expression for electron-hole friction force 
one should linearize the collision integral 

4 J M/(g){/e(p-q)/h(pi+q)[l-/e(p)][l-A(pi)]- 



7e (p)/h (Pi ) [1 - fe (P - q)] [1 - A (Pi + q)] } - 



(Al) 



where the non-equilibrium distribution functions for elec- 
trons and holes are determined via Eqs. © and g]), the 
factor of 4 arises from 2 possible spins and valleys for 
particle pi. W{q) is the Coulomb scattering probability 
which can be written down as 



Wiq) 



27r 



K{q + qTF) 
X 5{p + pi 



\{Up\Up')\^\{Up,\Up,^)\^ X 

|p-q| + |pi+q|) (A2) 



where qxF = 4aT/t>Fln(l -I- e^=/^)(l 



IS 



the Thomas-Fermi momentum describing screening in 



grapheni 



.24.34 



|(Mp|Mp' 



(1 + COS 9 pp') 12 is the ma- 



trix element of two Bloch functions in honeycomb lattice. 
In Eq. (jA2p the Thomas-Fermi momentum could be re- 
placed by the reciprocal gate-dielectric thickness for the 
qualitative description of gate screening. 

In case of small velocities: p • Vg <C T, p • V/i ^ T the 
electron-hole collision integral can be transformed using 



the common linearization technique for Fermi-systems^: 

St{feJh} = 

4 j W'(g)/e(p)/h(pi)[l - /e(p - q)][l - /h(pi + q)]x 



q-(V/^-Ve)rf^qd^Pi 

T {2TThY 



, (A3) 



where the energy and momentum conservation laws were 
used. 

Eq. (jASp . timed by q/2 and integrated over dFp is 
electron-hole friction force ^eh per unit area in the right 
hand side of the Euler equations ([SJ EI • After averaging 
over angle between q and V/i — Vg it is presented as 



where 



eh = I T^(g)/e(p)/h(Pl)[l-/e(p-q)][l-A(Pl+q)]x 

4g2 (fqd'^pid'^p 



T {2Trhf 



(A4) 



Ii>ii^ it was shown that almost collinear vectors p, 
Pi and q give the leading contribution to the collision 
integral. This facts originates from the linearity of the 
carrier spectrum, i.e. particles moving with the same 
velocity and the same direction interact infinitely long. 
Following the technique, described in^^, the quantity f3eh 
in the collinear limit is presented as 



rVhi(l/a) 



nOO pOD pP 

/ dP / dPi / F{P - z,)F{Pi + Zh) [l-F{P-Q- ze)] 



[1 - F(Pi +Q + zh)] y/PPi{P-Q){Pi+Q) 



Q^dQ 



i\Q\-4a\nFize)F{-Zh)y 



r, (A5) 



where F{x) = (l-fe^)"""^, z^ = fJ-e/T, Zh = pihlF , plasma. The momentum relaxation rate is given by: 
and the term Inl/a originates from the corrections to 
the electronic sperctrum due to electron-electron inter- 
actions. In this form the friction coefficient f^eh can be 
simply computed numerically. 



--/= / S,M/(g)(l-cosV)l("pl"p')l 



2 JV_ 

{2TiKf 



(A6) 



B. Friction caused by charged impurities and 
phonons 



Calculation of friction force caused by charged impuri- 
ties can be drawn analytically in the limit of monopolar 



where Si is the sheet density of charged impurities. The 
Coulomb scattering probability W{q) is given by eq. (jA2[) 
with the following delta-function: S{pvF — p'vp)- Here 
we put into consideration only self-screening caused by 
carriers in graphene. 

To calculate the friction force per unit area fg,; one 
should integrate the momentum relaxation rate (|A6[) 
timed by p with the non-equilibrium part of distribution 
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function Sfe 



de 



fe. = J P^P ■YeTpjdTp. (A7) 



In monopolar case one can set dfe^o/ds = —S{pvF — Me) 
and perform trivial integration: 

7re*SeSiVe f^^ SIT? OppidOppi 



nvj,K^ Jo (2 sin(6lpp72) + 4q)2 

where the dimensionless integral actually depends only 
on the permittivity k of the environment. 



We can use the same technique to derive the expression 
for the friction force, caused by phonon scattering. We 
start from the expression 



p^ph 



D^Tp 



derived in^, and, using Eq. (|A7|) . arrive at the following 
friction term: 



^eph 



(A8) 
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